Nonequilibrium Dynamic Phase Diagram for Vortex Lattices 
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The new dynamic phase diagram for driven vortices with varying lattice softness we present here 
indicates that, at high driving currents, at least two distinct dynamic phases of flux flow appear 
depending on the vortex-vortex interaction strength. When the flux lattice is soft, the vortices flow 
in independently moving channels with smectic structure. For stiff flux lattices, adjacent channels 
become locked together, producing crystalline-like order in a coupled channel phase. At the crossover 
lattice softness between these phases, the system produces a maximum amount of voltage noise. Our 
results relate spatial order with transport and are in agreement with experiments. 

PACS numbers: 74.60. Ge 
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Nonequilibrium problems involving elastic lattices and 
disordered media, such as the nature of depinning tran- 
sitions or the behavior of rapidly driven lattices, appear 
in a wide variety of systems including superconducting 
vortex lattices, charge-density waves, and solid-on-solid 
friction. Much recent interest has been devoted to the 
motion of a vortex lattice (VL) across a disordered sub- 
strate under applied currents both at and well beyond 
depinning. The reordering of a rapidly-driven VL is sup- 
ported by simulations as well as neutron scattering 
and decoration experiments , but the nature of the 
order that appears remains a subject of debate. In addi- 
tion, the relationship of the voltage noise observed near 
depinning to the microscopic vortex motion at higher cur- 
rents has not been addressed. 

Recent work on the behavior of a VL driven by a large 
current produced conflicting pictures of the VL order, 
ranging from a crystalline structure Q to a moving smec- 
tic state H . In Ref. ||, the VL does not fully recrystal- 
lize, but instead enters a moving Bragg glass state with 
channels. Moreover, Ref. || describes the VL in terms 
of independently moving channels of vortices with over- 
all smectic order. Other theories also focus on channels 
of vortices jl(|. Smectic structure factors (5(q)) of the 
VL were observed in simulations of vortices moving over 
strong pinning ||,^| , in agreement with Ref. || . 

Very recent decoration experiments [fl2f produced both 
crystalline-like and smectic order of the moving VL, de- 
pending on the magnitude of the applied magnetic field. 
Smectic order appears at low fields, when the vortices in- 
teract weakly, and crystalline-like peaks in S(q) appear 
at high fields, when the vortex interactions are stronger. 
This suggests that the softness of the VL is important in 
determining the vortex behavior at high driving currents. 

In this paper, we propose a new dynamic phase di- 
agram in which both smectic and crystalline-like order 
appear as the VL softness is varied. Using simulations 
of current-driven vortices, we clearly define regions of the 
phase diagram based on <S(q), V(T) curves, voltage noise, 
velocity distributions, direct observation of the lattice, 
and defect density calculations. We compute experimen- 



tally relevant voltage noise spectra at all currents 
from depinning to high drives, and find evidence for a 
washboard frequency in the stiff VL as well as broad- 
band noise in the plastic flow phase. As the VL softens, 
we observe a novel crossover to a regime in which smec- 
tic order is never destroyed even at high drives. At the 
crossover, the amount of voltage noise generated by the 
system during depinning is maximum, indicating that ex- 
perimental voltage measurements can probe the VL or- 
der. Our results are in excellent agreement with experi- 
ments, and have implications for noise measurements in 
the peak effect regime |l| . 

We model a transverse two-dimensional slice (in the 
x—y plane) of a T = superconducting infinite slab con- 
taining rigid vortices that are parallel to the sample edge 
(H = Hz). The vortex- vortex repulsion is correctly rep- 
resented by a modified Bessel function, K\(r/X), where A 
is the penetration depth. The vortices are driven through 
a sample with periodic boundary conditions, filled with 
randomly placed columnar defects, by a uniform Lorentz 
force fdX, representing an applied current. The colum- 
nar pins are non-overlapping, short-range parabolic wells 
of radius £ p = 0.30 A. The maximum pinning force, f p , 
of wells has a Gaussian distribution with a mean value 
of f p — 1.5/o and a standard deviation of 0.1/o, where 
fo = *§/87r 2 A 3 . The pin density n p = 1.0/A 2 is higher 
than the vortex density n v = 0.7$o/A 2 . We concentrate 
on samples 36A x 36A in size, containing 864 vortices and 
1295 pins. We check for finite size effects using samples 
that range in size from 18A x 18A to 72A x 72A and contain 
up to 2484 vortices and 3600 pins. 

The overdamped equation of vortex motion is fj = 
f™ + + id = T]Vi , where the total force fi 
on vortex i (due to other vortices f™, pinning sites 
f- , and the driving current f^) is given by f; = 



id- Here, O is the 
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Heaviside step function, (vj) is the location (velocity) 
of the zth vortex, rj^ is the location of the fcth pinning 
site, £ p is the pinning site radius, N p (N v ) is the num- 
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ber of pinning sites (vortices), = (r, — r ^ ) / 1 — Tj\, 

f ife = ( ri - r W)/| ri -rW|. 

Using the monotonic relationship between the shear 
modulus C66 and the dimensionless prefactor A v of the 
vortex- vortex interaction term, we model VLs with vary- 
ing degrees of softness by changing A v over 3.5 orders of 
magnitude, from A v = 0.001 to A v = 6.0. This is in con- 
trast to previous simulations ||-||,[lj,|l5| that considered 
only an extremely soft VL. For each value of A Vl we sim- 
ulate a voltage-current V{I) curve by initially placing the 
vortices in a perfect lattice, slowly increasing the driving 
force Id from zero to 3.0/o and measuring the resulting 
voltage signal V x — 1£ /N v . 

In order to identify the phase boundaries, each time 
the driving current id increases by O.O8/0, we halt the 
increases in f^, check that the voltage signal V x is sta- 
tionary over time, and then collect detailed velocity and 
position information for a long time interval, 2 x 10 s 
molecular dynamics (MD) steps, at a single current. For 
a 36A x 36A sample containing 864 vortices, each V(I) 
curve requires 8 x 10 6 to 10 7 MD steps to complete, cor- 
responding to about 10 days of CPU time on a SPARC 
Ultra. At each current, we are able to compute accurate 
voltage noise spectra S{uj) with a frequency window ex- 
tending into relatively low frequencies on the order of the 
vortex time-of-fiight across the sample. We also deter- 
mine the voltage noise power Sq in one frequency octave, 
So = QduS(uj), where S(u) = J V^e- 12 ^ dt\ 2 . 
The units of frequency are inverse MD time steps. Here, 
ujq = 0.027 and lj% — 0.054 were chosen to fall in the 
middle of our frequency window. The results are not 
affected if nearby values are used. To measure the or- 
der in the VL, we use the Delaunay triangulation of the 
instantaneous vortex positions to determine the fraction 
of six- fold coordinated vortices, Pq, and we compute the 
structure factor, 5(q) = | J2?=i e l ^\ 2 /N v , of the VL. 

The depinning of the VL appears in our simulated V(I) 
curves, shown in the inset to Fig. 0(a). Since a soft VL is 
able to deform and take better advantage of the available 
pinning sites, as A v decreases the depinning transition 
shifts to higher driving forces f^, from f^ = 0.48/o at 
A v = 6.0 to i d = 1.36/o at A v = 0.01. The depinning 
also becomes more abrupt for softer VLs, producing a 
peak in dV/dl that grows in magnitude with lower A Vl 
as seen in Fig. 0(a). The fact that a soft VL has a higher 
depinning current agrees well with experiments p2| , ^3[ . 

The VL also produces the greatest voltage noise at 
currents just above depinning, when the vortices are in a 
plastic flow regime. As seen in Fig. 0(b), for each value 
of A v the noise power Sq reaches its peak value S max at a 
driving current Lj above the current at which the peak in 
dV/dl occurs. The current id at which the peak occurs is 
independent of sample size. The spectra S(lo) near S'max 
are of the form S(lu) ~ as illustrated in Fig. 0(A). 
As A v decreases and the VL softens, we observe a peak 
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Driving Force f d /f 
FIG. 1. (a): dV/dl curves for vortex-vortex interactions 
(right to left) A v = 0.01, 0.1, 0.25, 0.5, 0.75, 1, 1.5, 2, 3, 4, 5, 
6. The peak in dV/dl increases in magnitude with decreasing 
Av Inset to (a): Voltage-current V(I) curves for the A v 
values listed above. The depinning transition shifts to higher 
driving currents f^ and becomes more abrupt as A v decreases. 

(b) : Voltage noise power So versus id for the A v values listed 
above. In each case the noise power peaks during the plastic 
flow regime. Inset to (b): Maximum noise power S miH as 
a function of A v . The peak value of S max corresponds to 
A v « 0.75, the same A v at which the system crosses between 
smectic and coupled channel behavior at high drives, (c): 
Fraction of six-fold coordinated vortices Pg as a function of 
id for values of A v listed above. At zero drive and strong VL 
coupling, Pa ~ 1 since the VL is field-cooled. The lowest value 
of Ps at each A v corresponds exactly to the peak in dV/dL 
The VL eventually reorders to Pe ~ 1 only when A v > 1. 
For A v < 0.75, P$ saturates at a value below 1. Inset to 

(c) : Velocity distribution functions D(v x ) for (1) plastic, (2) 
smectic, and (3) coupled channel phases. 

in S max at A v ~ 0.75, shown in the inset to Fig. 0(b). It 
is important to point out that this result agrees well with 
experiments conducted in the peak effect regime ]l3]| , in 
which an observed peak in voltage noise power 
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FIG. 2. Structure factor S(q) (a,c,e) and Delaunay trian- 
gulations (b,d,f) of the VL for: (a,b) A v = 1.5, i d = 1.12/ , 
in the plastic flow regime. S(q) is liquid-like, and the VL 
is filled with defects (marked with circles). (c,d) A v — 0.5, 
id = 3.04/q, in the uncoupled channel regime. S(q) has smec- 
tic peaks, and the defects in the VL are oriented transverse 
to the ^-direction driving force. (e,f) A v = A, id = 3.04/o, 
in the coupled channel regime. 5(q) has slightly anisotropic 
crystalline- like peaks, and the VL contains almost no defects. 

with changing magnetic field is interpreted to occur due 
to the softening of the VL as the externally applied mag- 
netic field increases. 

We quantify changes in the VL structure by calculating 
the fraction of six- fold coordinated vortices P§, shown in 
Fig. |l|(c), and the structure factor S'(q), shown in Fig. [|. 
The vortices are initially field-cooled, so Pq « 1 at f^ = 0. 
At driving currents below the depinning transition, the 
VL relaxes into the pinning sites, causing a gradual de- 
crease in Pq. The relaxation is most pronounced for soft 
VLs with low values of A v . A plastic flow state appears 
just above depinning, producing a liquid-like structure 
factor 5*(q), shown in Fig. |2|(a). The Delaunay triangu- 
lation in Fig. |^(b) reveals the large number of defects in 
the VL in this regime. As shown in Fig. [I], we find that 
the peak in dV/dl corresponds exactly to the driving 
force at which Pq reaches its lowest value. Thus dV/dl 



can be used to probe the order in the VL. 

As seen by the rise of P e in Fig. |l](c) for A v > 0.1, 
the VL regains order as the driving force f^ is increased 
and the interactions with pinning sites become less im- 
portant. In stiff VLs with large values of A v , the defects 
quickly heal out and Pq returns to a value near 1 for 
driving forces f^ not much larger than the pinning force 
f p . As A v is lowered and the VL softens, higher and 
higher drives f^ must be applied to bring P@ back to 1. 
In Fig. 0(c), for A v < 0.75, we see that Pq saturates at a 
value less than 1, and does not increase further even after 
applying driving forces id considerably larger than those 
shown in the figure (f^ ~ 6/0) • The permanent presence 
of a significant number of defects represents an impor- 
tant change caused by the VL softness in the behavior of 
the system at high driving currents. 

Computations of the structure factor S(q), shown in 
Fig. ^, verify that two different types of VL order appear 
at high driving currents depending on the lattice softness. 
When A v is low and the VL is correspondingly soft, the 
structure factor S(q) [Fig. ||(c)] has a smectic appear- 
ance, with well-defined peaks only along the q x = axis. 
A Delaunay triangulation of the moving VL, shown in 
Fig. H(d), reveals that the vortices are flowing in chan- 
nels oriented in the a;-direction (parallel to f<j) and ap- 
proximately regularly spaced in the y-direction. (This is 
in contrast to the plastic flow regime, in which vortices 
that remain pinned act as barriers and result in vortex 
motion transverse to the direction of flow, preventing the 
formation of longitudinal, straight channels.) Vortices in 
adjacent channels interact so weakly that each channel 
is able to move independently of nearby channels. All 
of the defects remaining in the VL have their Burgers 
vectors aligned transverse to f^, and each channel passes 
through at least one defect. Thus the channels can easily 
slip past each other as the defects glide, resulting in an 
uncoupled channel phase. 

If A v is large and the VL is correspondingly stiff, S(q) 
has a crystalline-like structure, shown in Fig. ||(e), at high 
driving currents. The peaks in S'(q) are anisotropic, with 
slightly higher peaks at q x = than at q x ^ 0. For larger 
samples, the same behavior is observed in 5(q) except the 
peaks are better resolved. The VL is nearly defect-free, 
as seen in the Delaunay triangulation of Fig. ||(f). The 
vortices still move in channels oriented in the direction 
of drive, but these channels are now locked together by 
the strong vortex-vortex interactions. Thus, the system 
is in a coupled channel dynamic state. 

We summarize the transitions between different states 
of the moving VL when A v and f^ are varied in the phase 
diagram shown in Fig. || The boundary between the 
pinned and plastic flow phases reflects the increase in de- 
pinning current as A v decreases, noted earlier. The plas- 
tic flow regime, identified by its liquid-like S'(q) and by its 
bimodal vortex velocity distribution [inset to Fig. 0(c)], 
disappears above id ~ f p . For id > f P , the softness of 
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Driving Force fjf 
FIG. 3. Dynamic phase diagram for different driving cur- 
rents id and vortex- vortex interaction strengths A v . At very 
low drives, the VL is pinned. In the plastic flow phase, we ob- 
serve peaks in dV/dl (crosses) as well as peaks in the voltage 
noise power So (circles). When A v < 0.75, the VL flows in un- 
coupled channels for all high driving currents. For larger A v , 
the VL passes through a transition region (UC) in which the 
channels progressively couple, until reaching a reordered state 
at high drives. A: Typical voltage noise spectrum S(u) ~ 1/ui 
observed near S ma>; for A v = 0.1. B: Washboard frequency 
observed in the coupled channel regime for A v — 2. The mag- 
nitude of the narrow band signal decreases in samples larger 
than the 36A x 36A sample shown here. The letters A,B refer 
to inset plots; a,c refer to Fig. ^ plots. 

the VL determines the behavior of the system. When 
A v < 0.75, the structure factor has a smectic appearance 
for all drives id ^ fp, corresponding to an uncoupled 
channel (UC) flow regime. For 1 < A v < 3, the sys- 
tem enters this same smectic UC state immediately after 
leaving the plastic flow state. At slightly higher driving 
currents, however, the defect-filled VL enters a transition 
regime marked by the appearance of weak peaks in 5(q) 
at q x ^ 0. These peaks are significantly smaller than 
the smectic peaks at q x = 0, but are at least twice as 
large as the background observed in the smectic phase. 
Throughout the UC transition regime, the channels grad- 
ually couple as the number of defects in the VL drops. 
The unimodal vortex velocity distribution becomes less 
broad, as in the inset to Fig. 0(c). When A v > 4, the sys- 
tem goes directly into this transition state without ever 
passing through a truly smectic state. Eventually, for all 
VLs with A v > 1.0, the system reaches a coupled channel 
flow phase for large enough values of f^. The boundary of 
this phase is identified as the current at which Pq reaches 
a value of 1. At the transition A v ~ 0.75 below which the 
coupled channel phase never appears, the greatest max- 
imum noise power S'max is observed, indicating that the 
noise power could be used as an experimental probe of 



this transition. At A v ~ 0.75, vortex-vortex and vortex- 
pin interactions are nearly balanced, allowing the system 
to sample the largest number of metastable states. The 
boundaries shown on the phase diagram mark crossover 
points rather than sharp phase transitions. 

In the coupled channel regime, we find a washboard 
frequency |l6| in the voltage noise signal, shown in 
Fig. ||(B), corresponding to the time required for a vor- 
tex to move one lattice constant. The magnitude of this 
washboard signal decreases when the system size is in- 
creased, indicating that the signal appears only when the 
region of the VL sampled is locked into a single domain. 
Thus experimentally it would be necessary to probe the 
voltage over a very small area of a sample to observe a 
washboard frequency, such as with local Hall probes. 

In conclusion, we have obtained a new vortex dynamic 
phase diagram as a function of lattice softness and driving 
current. At high driving currents two distinct phases of 
flux flow appear: in soft lattices, uncoupled channels with 
a smectic structure, and in stiff lattices, coupled channels 
with crystalline-like order. A signature of the crossover 
is observed in the voltage noise, which is largest at the 
transition between the two phases. In the coupled chan- 
nel phase, a washboard frequency appears in the voltage 
noise spectrum for small sample sizes. Our results are in 
agreement with experiments. 
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